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p ^ ■ Abstract. We study the linearized Fopl - von Karman theory of a long, thin 

rectangular elastic membrane that is bent through an angle 2a. We prove rigorous 
bounds for the minimum energy of this configuration in terms of the plate thickness 
a and the bending angle. We show that the minimum energy scales as er 5 / 3 a 7 / 3 . This 
scaling is in sharp contrast with previously obtained results for the linearized theory 
of thin sheets with isotropic compression boundary conditions, where the energy scales 
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1. Introduction 



Everyday experience tells us that thin elastic sheets crumple, when confined into a 
small volume, e.g. a sheet of paper "confined" by ones hands. Crumpling also plays an 
important role in the mechanical behavior of packaging material and in the dissipation 
of the energy of collisions by the "crumple zones" of automobiles. Crumpling is therefore 
a problem of much intrinsic interest, but understanding this behavior is complicated by 
the complex morphology of a typical crumpled sheet. 

Despite the complicated appearance of a crumpled sheet, the crumpling 
phenomenon is in itself very robust. It is easily observed in thin sheets made from 
a variety of materials, suggesting that it can be studied using simplified or idealized 
models that capture the essential features of thin elastic sheets. This approach leads 
one to consider a crumpled sheet as a minimum energy configuration for a simple elastic 
energy functional for thin sheets, viz. the Fopl - von Karman (FvK) energy. Using this 
approach of elastic energy minimization, the crumpling response is now understood as a 
result of the elastic energy of the sheet concentrating on a small subset of the entire sheet 
[| . The energy in a crumpled sheet is concentrated on a network of thin line-like 
creases (ridges) that meet in point-like vertices. Recent work has resulted in quantitative 



understanding of both the vertices |4], ||, ||, [7], || and the ridges || |10[. Scaling laws 
governing the behavior of crumpled sheets have been obtained in the physics literature 

§ 0, il- 

Minimum energy configurations for the FvK energy have also been studied in 
the context of the blistering problem, viz. the buckling of membranes as a result of 
isotropic compression along the boundary. The blistering problem is relevant to the 
delamination of thin films that are chemically deposited at high temperatures, as well 
as the mechanical behavior of micro-fabricated thin-film diaphragms [11 , |l^ . 

There is a considerable body of mathematical work focused on the blistering 
problem [|ll], [13], [14], [15|, [16], |17], HH. Upper and lower bounds have been obtained for 



approximations to the elastic energy [|TT], [13], |16[ , for the FvK energy [ 15|, [XT] and for full 
three dimensional nonlinear elasticity [P^] . The FvK energy and full three dimensional 
nonlinear elasticity give the same scaling for the upper and the lower bounds. This 
yields a rigorous scaling law for the energy of a blister. As we discuss below, the scaling 
laws for the blistering problem are different from those for the crumpling problem, even 
for the same energy functional. This indicates that the energy minimizing configurations 
for the FvK energy show an interesting dependence on the boundary conditions. 

The minimization of the FvK elastic energy is an example of a non-convex 
variational problem that is regularized by a singular perturbation []I9|, pD . It is well 



known that this can lead to a variety of multiple-scale behaviors including energy 
concentration and/or small scale oscillations in the minimizers |[21|| . 

Multiple scale behaviors, both micro structure and singularities, are ubiquitous in 
condensed matter systems. The crumpling phenomenon appears to be a particularly 
simple and tractable instance of multiple scale behavior. For this reason, there has been 
much recent interest in physics literature about the nature of the crumpling phenomenon 
22) , p3j [p. ||; |24| . Asymptotic analysis |J and scaling arguments || show that the ridge 



energy scales as 

£ r ~ a^L 1 / 3 , 

where a is the thickness of the sheet, and L is the length of the ridge. Ben-Amar et al 
]]J and Mahadevan et al ||] show that the energy of a <i-cone scales as 

£ c ~ a 2 log(L/a), 

where a is the radius of the core associated with the vertex. It is clear that the ridge 
energy is asymptotically larger than the vertex energy as a — > 0. However, which 
of these two energies is important for a given sheet depends on the relation between 
the nondimensional thickness e = a/L of the sheet, and the "crossover thickness" 
e* = a/L which is determined by setting £ r = £ c . e* in turn depends on the values 
of the multiplicative constants for these scaling laws. 

These multiplicative constants cannot be determined by scaling arguments. 
Mahadevan et al estimate the constant for the vertex energy numerically using an 
ansatz for the shape of a sheet near the vertex, and they find that 

£ c w 1000V 2 log(L/a), 
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where is the complement of the tip angle of the cone || Based on Lobkovsky's 
work |J, Boudaoud et al |7J estimate the value of the constant in the ridge energy as 1. 
This implies, for the experiments in Refs. || and O, the vertex energy dominates the 
ridge energy. 

It would be useful to prove these scaling laws, and determine the multiplicative 
constants rigorously, that is, in an ansatz-free manner. In this paper, we propose a 
model problem that yields structures analogous to a single ridge in a crumpled sheet. 
We prove a rigorous lower bound (with a numerical value for the multiplicative constant) 
for the elastic energy in our model problem. This is a step toward rigorously proving 
the scaling law of Lobkovsky et al for the ridge energy ||. We also discuss how the 
techniques in this paper can be extended to prove similar results for a "real" crumpled 
sheet, as opposed to our model problem. 

This paper is organized as follows - In Sec. |2j, we describe the problem of interest, set 
up the relevant energy functional and determine the appropriate boundary conditions. 
We also rescale the various quantities to a form that is suitable for further analysis, and 
recast the problem in terms of the rescaled quantities. In Sec. ^ we prove our main 
result, viz. a lower bound for the elastic energy for our boundary conditions. We present 
a concluding discussion in Sec. f|. 

2. The Elastic energy 

We are interested in a minimal ridge, i.e., the single crease that is formed when a long 
rectangular elastic strip is bent through an angle by clamping the lateral boundaries to 
a bent frame. This situation is depicted in Figure [I|. 

From the symmetry of the problem, it is clear that we only need to consider one half 
of the strip. This is represented schematically in Figure We will use (the material) 
coordinates (x, y) on the reference half strip \x\ < L, < y < oo. Also, we associate a 
{u, v, w) coordinate systems in space, so that as y — > oo, the half-strip is asymptotically 
in the w = plane, as depicted in Figure |2|. The grid in the figure is generated by 
the lines x = constant and y = constant that are straight in the reference (material) 
coordinates, w represents the out of plane deformation, and the in-plane coordinate 
are chosen so that the u and the v axes are asymptotically parallel to the x and y axes 
respectively as y — > oo. Since the sheet is bent through an angle 2a, as y — > — oo, the 
sheet will lie in the plane w = v tan(2a). The symmetry of the two halves implies that 
the the line y = maps into the plane w tan a + v = 0, which bisects the angle between 
the planes w = (the asymptote as y — ► oo) and w = v tan(2a) (the asymptote as 
y — > — oo). The symmetry between the two halves also implies that the tangent to the 
lines x = constant at y = should be perpendicular to the plane w tana + v = 0. 
Consequently u y = and w y = v y tan a at y — 0. 

A mathematically justified way to obtain the elastic energy of the deformed sheet is 
to treat the sheet as a three dimensional (albeit thin) object and use a full nonlinear three 
dimensional elastic energy functional for the energy density. This approach however does 
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Figure 1. A minimal ridge. The boundary conditions are given by a frame (the 
thick solid lines) bent through an angle. The sheet is essentially flat outside the region 
bounded by the two dashed curves, and the bulk of the energy is concentrated in this 
region. 



A w 




v = v Q (x) 



Figure 2. A schematic representation of our coordinate system and the boundary 
conditions imposed on the sheet. The grid is given by the lines x — constant and 
y = constant. The two dashed lines are the tangent to the curve x = and the line 
w = wq{Q),u = respectively, and the angle between theses lines is a. 
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not take advantage of the "thinness" of the sheet. In particular, we would like to treat 
the thin sheet as a two dimensional object. This will greatly reduce the complexity of 
the problem. The derivation of reduced dimensional descriptions of thin sheets has a 
long history. There is a classical theory for thin elastic sheets built on the work of Euler, 



Cauchy, Kirchoff, Fopl and Von Karman 25, 26 



Many modern authors have considered the problem of deriving a reduced 
dimensional theory pTj from three dimensional elasticity in a mathematically rigorous 
fashion through T - convergence [[28], |29| . T - convergence is the appropriate notion for 
convergence in variational problems. In the context of thin sheets, roughly speaking, 
finding the T - limit amounts to identifying an appropriate two dimensional energy 
functional whose minimizers give the limiting behavior of the minimizers of the full 
three dimensional energy in the limit the thickness o — > 0. This problem hasn't been 
solved in its entirety. 

There are two situations for which reduced dimensional theories have been derived 
rigorously as T - limits of full three dimensional elasticity. Membrane theories |3(], 31], ^ 
are applicable in situations where the stretching is essentially uniform through the 



thickness of the sheet and bending theories [|33|] are applicable in situations where the 
the sheet is essentially unstretched, and all the elastic energy is due to strains that are 
first order in the thickness of the sheet. 

Neither of these theories are applicable for the minimal ridge. For the minimal ridge, 
both the stretching (membrane) energy and the bending energy are important. In fact, 
a ridge is a result of the competition between these two energies. So we turn to the 
classical Fopl - von Karman ansatz for an appropriate definition of the elastic energy 
f25| , [26| . Although the FvK energy is not rigorously derived from three dimensional 
elasticity, it can be thought of as the sum of the membrane and the bending energies 
that have been derived rigorously in different limiting situations. We will further assume 
that, for the minimal ridge, the deflections \w\ ) \u — x\ and \v — y\ are small compared to 
the natural length scale L. The strains are of the order of the square of the maximum 
deflection divided by L and they are small. After some rescaling, the energy of the 
deformed sheet is given by the linearized FvK energy 

1 



1 



(U X +W X -1) + ~(V X + Uy + 2W X Wy) +(Vy+W y — iy 



dxdy 



w 



dxdy. 



(1) 



where x and y are reference coordinates on the sheet, u and v are in-plane coordinates, 
w is the out of plane displacement and a is the thickness of the sheet. The integrand 
for the first integral is given by the squares of the linearized strains, 



U x -\- W x 1, Ixy OVr 



W X Wy + ~(V X + Uy), =V V + W y -1. 



The blistering of thin films is also described by the elastic energy in (fj). A similar 
energy also describes multiple scale buckling in free elastic sheets (i.e. sheets that are 
not forced through the boundary conditions) that are not intrinsically flat |34 . 
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The difference between the blistering problem and a minimal ridge is in the 
boundary conditions, which we describe below. Since the strains are assumed to be 
small, v y + Wy ~ 1. If the bending half-angle a < 1, so that tana « a, we get 
w y = a at y = 0. Since the deformation goes to zero far away from the bend, we have 

| it — x\ — > 0, \v — y\ — > 0, w — > as y — > oo, 

Also, the lateral boundaries at x — ±L are clamped to the frame. Therefore, 

u = x, v = w = 0, at x = ±L. 

Note that the appropriate boundary condition for the minimal ridge at y = is a free 
boundary condition, subject to the symmetry requirement w + v tana = 0. We are 
going to replace this free boundary condition with a Dirichlet boundary condition 

v = v (x),w = w (x) at y = 0, 

where we will leave the functions Vq(x) and Wq(x) unspecified, except for a size condition. 
Defining 

a = maxO (x)|, \w (x)\), 

[—L,L] 

we impose the size condition by requiring that a be "small" . 



2.1. Rescalings 

The relevant small parameter in the problem is the non-dimensional thickness of the 
sheet, e — a/L. Following Lobkovsky Q, we introduce the rescaled coordinates and 
displacements by 

x = LX, y = a 1/3 L 2/3 Y, 

and 

w = a^L^W, v = y + a^L^V, u = x + a^L^U. 

Since a,L,x,y,u,v,w all have dimensions of a length, it is clear that the rescaled 
quantities X, Y, U, V, W are all dimensionless. With these rescalings, the dimensionless 
energy I = <r~ 5//3 L~ 1//3 Z is given by 



I(U,V,W) 



-2/3 



(U x + W 2 X ) 2 + -— (V x + U Y + 2W X W Y ) 2 + e~ 4/3 (W + W, 



2\2 
YJ 



2 

+ [W 2 Y + 2e 2/3 W 2 Y + e i/z W 2 xx ] dXdY. (2) 

Our quest for rigorous scaling results for the energy X reduces to the following - Show 
that the rescaled energy J, of a minimizer (U*, V*, W*), is bounded above and below by 
positive constants uniform in the dimensionless thickness parameter e, as e — > 0. 
Setting 

W (X) = a-^L-^woiLX), V (X) = a^L^VLX), 

6 



the rescaled quantities satisfy the boundary conditions 

V = V (X),W = W (X),W Y = a atr = 0, 

and 

U -> 0,V -> 0, W -> asF^oo, 
We have the lateral boundary conditions 

£7 = y = W = at X = ±1. 

Also, the deformation at Y = satisfies |Vbpf)| < A and |Wq(X)| < A where 
A = a-^L-Wa. 

3. Lower Bound 

In this section, we prove a lower bound for the linearized Elastic energy in Eq. flU), 
provided that the length scale a associated with the boundary conditions satisfies a size 
condition. 

Theorem 3.1. T(u,v,w) is as defined in Eq. ^). There is a constant b > such that, 
for all u G H l ,v G H 1 and w G H 2 H W 1 ' 4 satisfying the boundary conditions 

\u — x\ ^ 0,\v — y\ ^ 0,w ^ as y — > oo, 
u = x, v = w = 0, at x = ±L, 
v = v (x), w = w (x) at y = 0, 

and the size condition 

maxOoO)|, \w {x)\) = a< ba 1/3 L 2/3 a 2/3 , 

[-L,L] 

we have the lower bound 

I(n,ii,«;)>^ T /W. 
5 

Remark. We will prove the theorem by proving the scaled version of the statement, 
viz., the rescaled boundary conditions along with the rescaled size condition 

max(|K,(X)|, \W (X)\) = A < ba 2 ' 3 , 

imply that 

T 2a'/ 3 

I > . 

~ 5 

Remark. The hypothesis for the theorem includes a size condition on the displacement 
at y = 0. This is somewhat unsatisfying, as it is not a priori obvious that a ridge in a 
"real" crumpled sheet will satisfy this condition. 
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In our search for a lower bound, we can assume I(U, V, W) < oo w.l.o.g. From 
Eq. (Q), it follows that for a > 0, W G H 2 . By the standard trace theorems, the 
boundary conditions W(X, 0) = Wq(X) and Wy(X,0) = a are therefore assumed 
pointwise for almost every X. 

Since U = at X = ±1, it follows that /i^tf* + W^dX = Wj-dX. Using this 
together with Jensen's inequality we get 



I(U, V,W)> J [(U x + Wlf + W YY ] dXdY 



> 



i \ 2 /• j 

W x dX) + 
i / •/ -i 



(3) 



Thus the functional 

E(W) 



i \ 2 /-J 

-l / J-i 



dy 



bounds I(U, V, W) from below. Hence it suffices to obtain a lower bound for E with the 
given hypothesis to prove the theorem. 
Let Eb and E s denote the quantities 



E b = / WyydXdY, 



and 



E s 



W x dX | dY. 



"OO ^ 

'o 2 

Although i£f, and are only lower bounds for the "true" bending and the stretching 
energies I b and I s , we will nonetheless call Eb and E s the bending and the stretching 
energy for convenience. 

For every X, we define 



and for every Y, we define 



t(Y) 



Wyy{X,Y)dY 



W x dX. 



-i 



Since E < oo, p(X) > (a.e.) and r(Y) < oo (a.e.). 

For any function / depending on X and Y, let ||/||y denote f_ ± \ f(X,Y)\ 2 dX, so 

that, 

r(Y) = \\W x \\ 2 y. 

p(X) is a "local" (in X) measure of the bending energy, and [p(X)] _1 can be thought of 
as the bending energy density in X that is obtained by integrating out the Y dependence. 
We can also think of p(X) as the length scale associated with the bending energy as a 
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function of X, viz., we expect that the bending energy density in Y decays rapidly for 
Y/p(X) ^> 1 (See Fig. [I]). Likewise, [r(F)] 2 is a local (in Y) measure of the stretching 
energy. In terms of p(X) and t(Y), the bending an the stretching energies are given by 



We now begin our proof of the theorem. The idea behind the proof is to show that 
the stretching energy E s can be bounded from below by a negative power of the bending 
energy E b , so that the total energy E s + E b tends to +00 as E b — > and E b — > 00. This 
ensures the existence of a positive lower bound for E (and consequently also for J). 

Lemma 3.2. For every Y, we have the inequality 



Proof. If p(X) > 0, W(X, .) is a C 1 function by the Sobolev Embedding theorem and the 
boundary conditions imply that W(X, 0) = Wo(X) and Wy(X, 0) = a. Consequently, 
for such an X, 



We will estimate W 2 (X, Y) from this equation using the elementary inequalities 




\\W\\ Y > a 2 Y 2 - 2Y 3 E b - AA 2 . 




a-b\ 2 > (1 



S)\a 



1-5 

5 



M 2 



and 




W (X)\ 2 < A 2 . 



By Jensen's inequality 





Y 3 



Combining these estimates, we get 



W 2 (X,Y) > (1 -5i)a 2 F 2 



1 - 



(1 + 62) 



Y 3 



l + S 2 

$2 



A 2 



5 



+ 
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for all positive 5i and 8%. In particular, setting S\ = 1/2 and 82 = 1 yields 

W 2 {X ) Y)>\a 2 Y 2 -^- ) -2A 2 . (4) 

Integrating this inequality in X we obtain 

\\W\\ 2 Y >a 2 Y 2 -2Y 3 E b -AA 2 , (5) 
proving the lemma. □ 

Our proof is based on demonstrating that a small bending energy Eb will lead to 
a large stretching energy This idea is quantified by the following lemma where we use 
Eq. (^) to estimate the stretching energy E s from below in terms of the bending energy 
E h . 

Lemma 3.3. Let 

/ 4E b A x 2 

There is a constant fi* > such that, if /j < the stretching energy satisfies 

a 14 

K. > 



5 • 729 • El ' 



Proof. We have, 



7T _ 



r (y) = > — ||W||^ 



by the Poincare Inequality. Eq. (^|) now yields 



r(Y)>^[a 2 Y 2 -2Y 3 E b -4A 2 }. 

We will now extract the appropriate scaling dependence of t(Y) and E s on a and Eb- 
Setting a 2 Y 2 = 2Y 3 Eb, we deduce that a characteristic scale Y for Y is given by 



Y 



a 



2E h 



Rescaling Y in terms of Y, we obtain 



t(Y)> 



a Ti 
16E 2 



Y 
Y 



where /x is as defined above, i.e. 



AE h A 



a' 
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Making the change of variables Y = zY, and using r(Y) > 0, we see that 

7T 4 a 14 iT(/i) 



E s > J -[r(Y)fdY> 



512E, 



5 ' 



where 



/ [max(z 2 (l — z) — fi, 0)] 2 dz. 
Jo 
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K(fi) is clearly a continuous function of fi and -ft'(O) = 1/105. Since 7r 4 i^(0)/512 > 
1/(5 ■ 729), there is an //* > such that for all // < /i*, rr 4 K(fi)/512 > 1/(5 ■ 729). The 
lemma follows. □ 

We can now prove the theorem. 

Proof. Let b = f The hypothesis imply 

2/3 



If -Efe > 2a 7 / 3 /5, there is nothing to prove. Therefore, we can assume that E b < 2a 7/3 /5. 
Then, it follows that 

Consequently > a 14 / (5 ■ 729 ■ El). Minimizing E s + E b , we see that 

a 14 „ . 2a 7 / 3 



/ > E s + E b > r + 

- s - 5 ■ 729 ■ El ~ 

The theorem follows by "undoing" our rescaling to express the results in terms of 
u, v, w, a and X. □ 

Remark. In the appendix, we show that 

- 105 6 

Therefore, we can choose \i* = 0.048 and b = 0.13. We see that b is not exceedingly 
small. Rather it is 0(1). Also, it is not the best constant for this theorem, and we can 
get a better value by optimizing our choices for the constants (Si and 82 which we set 
to be 1/2 and 1 respectively). 



4. Discussion 



We have proved a rigorous lower bound for the elastic energy of a ridge that is 
consistent with the results obtained by Lobkovsky et al [0 using scaling arguments, 
and by Lobkovsky using matched asymptotics. In order to prove these scaling laws 
rigorously, we also need analogous upper bounds that are consistent with the same 
scaling. This approach has been used successfully for a variety of other variational 
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problems |35], [T|, |T7|]. The upper bound is usually obtained by constructing a test 
solution that yields the bound. One is often guided by scaling arguments in the 
construction of the appropriate test solution. This is in contrast to the lower bounds, 
where one needs to obtain functional analysis type inequalities that captures the 
competition between distinct energies in the problem (e.g. lemma |3.3| ). It is this 
competition that determines the scaling behavior of the problem. 

We have not constructed the upper bounds, because we believe that they will follow 
from the scaling ansatz in Lobkovsky's work |§. We also believe that the upper bound 
will scale in the same manner as the lower bound, thereby giving us a rigorous scaling 
law for the energy of a single ridge. 

As we remark before proving the theorem, our result is not directly applicable to 
the minimal ridge since we have the extra hypothesis 

We need to replace this restriction with a free boundary condition at y = subject to 
the symmetry requirement w tana + v = 0, and the boundary condition w y = v y tana. 

Despite this caveat, we claim that the analysis in this paper captures the 
essential features of the minimal ridge problem, namely the scaling in lemma |3.3| , and 
consequently the scaling law for the energy of the ridge. We will show that this is indeed 
the case by investigating the problem with the "true" boundary conditions elsewhere. 

A harder problem is to show that the scaling laws also hold for a real crumpled sheet, 
where the forcing is not through clamping the boundaries to a frame, but rather through 
the confinement in a small volume. In this case, there are interesting global geometric 
and topological considerations, some of which are explored in Refs. [[36] and | 57| . As 
Lobkovsky and Witten |TtJ argue, the boundary condition that the deformation goes 
to zero far away from the ridge implies that the ridges do not interact with each other 
significantly. The ridges can be considered the elementary excitations of a crumpled 
sheet. In our quest for upper and lower bounds, this translates to the statement that 
the test solutions for the upper bound can be constructed by piecing together local 
solutions near each ridge. Thus the upper bound is obtained relatively easily. The 
hard part is to show that the confinement actually leads to the formation of ridges, and 
that the competition between the bending and the stretching energy for this situation 
has the same form as in lemma |3.3| . In this context, we expect that global topological 
considerations, as well as the non self-intersection of the sheet will play a key role in the 
analysis, as they do in the analysis of elastic rods (one dimensional objects) ||38|| . 

As we note above, the blistering problem is described by the same elastic energy 
(Eq. (fj])), but with different boundary conditions. Our results show an interesting 
contrast with results for the blistering problem. Ben Belgacem et al. have shown that 
fl"7| , for an isotropically compressed thin film, the energy of the minimizer satisfies 

c\ 3/2 aL <1< CX 3/2 orL, 

where L is a typical length scale of the domain, and A is the compression factor. A 
construction for the upper bound strongly suggests that the minimizers develop an 
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infinitely branched network with oscillations on increasingly finer scales as cr — > 0. In 
contrast, our results indicate that the energy of a minimal ridge satisfies 

ca^L 1 ' 3 < J < Ca^L 1 ' 3 , 

and the energy concentrates in a region of width cr 1//3 L 2//3 . This shows that the nature of 
the solution of the variational problem for the elastic energy in ([I]) depends very strongly 
on the boundary conditions. In particular the very nature of the energy minimizers is 
different for the two problems - For the blistering problems, as o — > the minimizers 
develop a branched network of folds refining towards the boundary. The minimizers 
therefore display the problem of small scale oscillations |39|, |2l| . The minimal ridge 
problem on the other hand displays the concentration phenomenon [^, |2l| as a — > 0, 



with the energy concentrating on a region of width ~ a x l 3 L?l 3 . It would be interesting 
to explore this issue further, and in particular, to understand the mechanisms by which 
the boundary conditions determine the nature of the minimizer. 

Appendix 

In this appendix, we prove the inequality 



A 

Jo 



K(ji) = I \^{z\l-z)-^)] 2 dz>^-^. 

Proof. Let / > 0. Then, for / > /i, we have 

[max(/ - /i, 0)] 2 = f - 2nf + fi 2 > f - 2 M /. 
For < / < /i, we have 

[max(/- yU ,0)] 2 = 0>/ 2 -2/i/. 
Using these inequalities in the definition of K{p) we obtain 

K(fi) > / [z 2 (l -z)] 2 - 2fiz 2 (l - z)dz > — - 
Jo 
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□ 
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